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(54) Measuring absorption of a scattering medium 

(57) An apparatus for measuring absorption of a 
scattering medium (2), comprises a light source (10) for 
generating pulsed light, a light guide (12) for making the 
pulsed light incident on the scattering medium (2), a light 
guide (13) for receiving the light having propagated 
through the scattering medium (2), a signal detecting 
section (1 5) for temporally cutting a portion out of a light 
signal obtained by the light receiving section, a first 
arithmetic section (16) for acquiring measurement sig- 
nals respectively at plural different timings and for de- 
riving a time integration value and a mean optical path- 
length for each of the plural measurement signals ob- 
tained, and a second arithmetic section (17) for calcu- 
lating a change amount of the concentration of an ab- 
sorptive constituent. 
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Description 

BACKGROUND OF THE INVENTION 
5 Field of the Invention 

The present invention relates to a measuring method and apparatus capable of obtaining information concerning 
absorption inside a scattering medium like living tissues and the like, an Imaging method and apparatus, a fluoroscopic 
method and apparatus, a tomographic image taking method and apparatus, or a mammographic method and appara- 
10 tus. More particularly, the present invention concerns a method and apparatus for measuring a temporal change or a 
spatial distribution of concentration of a specific absorptive constituent in a specific portion inside a scattering medium 
and further concerns an absorption information measuring method and apparatus of scattering medium capable of 
non-invasively measuring a concentration distribution of a specific absorptive constituent inside a scattering medium 
having surfaces of various shapes. 
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Related Background Art 



There are four types of main methods for measuring or imaging a living tissue of a scattering medium by use of 
light, which are roughly classified into (a) one using translationally propagating light (ballistic photon), (b) one utilizing 
20 a coherent component by heterodyne detection thereof, (c) one utilizing scattered light, and (d) one utilizing a time- 
resolved gate. 

(a) The method utilizing the translationally propagating light (ballistic photon) is arranged to make pulsed light 
incident to a measured object and to extract a translationally transmitted component In output light by an ultra- 
25 high-speed gate. On this occasion, a flight time of the translationally transmitted light is shorter than those of 

scattered light and is shortest among the output light. Therefore, the translationally transmitted light can be detected 
by cutting out a light part first emerging from the measured object by the ultra-high-speed gate. This method has 
such an advantage that a high spatial resolution as in X-ray CT can be expected theoretically. However, almost 
all light is scattered multiply Inside the measured object being the scattering medium, and therefore, the transla- 
te tionally transmitted component exists little or, even if present, in a very small amount. As a result, the utilization 
factor of the light upon measurement becomes very poor and, considering the maximum light incidence intensity 
(approximately 2 mW/mm^) that can be injected into a living body, an extremely long time becomes necessary for 
imaging or measurement or the like, which is not practical. With a large object like the head or the arm the signal 
light is very weak and measurement is impossible in practice. 
35 (b) The method utilizing the coherent component Is arranged to split coherent high-frequency modulated light Into 
two components and make one incident to the measured object. Then the output light is made to interfere with the 
other high-frequency modulated light component to heterodyne-detect the coherent component included in the 
output light, and the coherent component is used for measurement. Specifically, since in this measuring system 
the flight time (or optical pathlength) and the flight direction of scattered light are disordered depending upon the 
40 degree of scattering, multiply scattered light does not interfere, and thus nearly translationally transmitted light is 
detected. Also in this case, therefore, the high spatial resolution is expected theoretically as in the previous ex- 
ample, but the utilization factor of light is very poor, so that it is very difficult to practically apply it to measurement 
or imaging or the like of living body. It is also practically Impossible to measure a targe object like the head or the arm. 
(c) The method for detecting the scattered light Is arranged to make continuous wave (CW) light, pulsed light, or 
modulated light Incident to the measured object and to detect the light (output light) having propagated as scattered 
In the measured object. On this occasion, where the pulsed light of a sufficiently shorter pulse width than the time 
response characteristic of the measured object is made incident thereto, it may be considered as incidence of 
Impulse light and the output light corresponding thereto is expressed by an impulse response. The measurement 
of this type is called as time-resolved spectroscopy (TRS). wherein internal information is calculated by measuring 
50 a time-resolved waveform of the impulse response. The method for making the continuous wave (CW) light incident 
Is called as CW spectroscopy and the measurement for making the modulated light incident as phase modulated 
spectroscopy (PMS). The former measures a time Integration value of the impulse response and the latter meas- 
ures a Fourier transform of the impulse response, i.e., a system function thereof. From the above discussion, it is 
understood that the time-resolved spectroscopy (TRS) obtains the signal containing more information than the 
other two methods. The phase modulated spectroscopy (PMS) for sweeping the modulation frequency and the 
time-resolved spectroscopy (TRS) of impulse response include nearly equivalent information. Since the continuous 
wave (CW) spectroscopy and the phase modulated spectroscopy (PMS), though greatly different in their signal 
processing methods, obtain signals Including substantially equivalent information, almost equal information can 
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be measured thereby. Since alt of these methods utilize almost all output light, their utilization factors of light are 
increased remarkably as compared with the foregoing two examples of (a) and (b). However, the scattered light 
spreads in the almost all area or in a considerably wide region inside the measured object, so that the above 
methods cannot be applied to quantitative measurements in a narrow, specific portion. In imaging or optical CT, 
s the spatial resolution is very poor and practical application is difficult. 

In contrast with the above, (d) the method utilizing the time-resolved gate uses a part of the impulse response 
waveform. For example, photons included in the initial portion of the time-resoived output signal (impulse response 
waveform) are those little scattered (photons little tarrying on the way), which are photons having propagated near the 

to straight optical path connecting a light incidence position (light incidence point) with a photodetection position (light 
receiving point). Accordingly, by utilizing the signal from the time-resolved gate, the photons having propagated near 
the above straight optical path and then detected can be extracted, thus improving the spacial resolution of imaging 
or the like. The extremity of this is the translationally propagating tight (ballistic photon) described previously. The time- 
resolved gate methods having been developed or proposed heretofore do not permit the quantitative measurements 

75 of an absorptive constituent inside the scattering medium. The reason is that no method is developed for analytically 
describing the signal obtained by the time-resolved gate method. Namely, there are no known methods for analyzing 
the signal obtained by the time-resolved gate method based on the photon diffusion theory and for utilizing it for quan- 
tification of absorptive constituent. Even if a method of analysis based on the photon diffusion theory is developed, 
there will remain the essential problem resulting from the basis placed on the photon diffusion theory; i.e., there will 

20 remain a problem that it is impossible to measure scattering media having various contours and those to which diffusion 
approximation cannot be applied. From the above, the time-resolved gate method can improve the spatial resolution, 
but has the serious problem that the internal information of an absorptive constituent or the like cannot be quantified. 

There are many references as to the time-resolved gate method as described above and typical examples thereof 
are listed below. It is, however, noted that none of them presents discussions or suggestions on the quantification of 

25 absorptive constituent as disclosed in the present invention. 

1 ) A. J. Jobtin, "Method of calculating the image resolution of a near-infrared time-of-flight tissue-imaging system. 
" Appl. Opt. Vol. 35. pp. 752-757 (1996). 

2) J. C. Hebden, D. J. Hall. M. Firbank and D. T Delpy. "Time-resolved optical imaging of a solid tissue-equivalent 
30 phantom," Appl. Opt. Vol. 34. pp. 8038-8047 (1995). 

3) G. Mitic, J. Kolzer, J. Otto. E. Plies, G. Solkner and W. Zinth, "Time-gated transillumination of biological tissues 
and tissue-like phantoms." Appl. Opt. Vol. 33, pp. 6699-6710 (1994). 

4) J. C. Hebden and D. T Delpy, "Enhanced time-resolved imaging with a diffusion model of photon transport," 
Opt. Lett. Vol. 19, pp. 311-313 (1994). 

35 5) J. A. Moon, R. Mahon, M. D. Duncan and J. Reintjes. "Resolution limits for imaging through turbid media with 

diffuse tight." Opt. Lett. Vol. 18, pp. 1591-1593 (1993). 

SUMMARY OF THE INVENTION 

40 An object of the present invention is to solve the problem of "incapability of quantitatively measuring the internal 

information of absorptive constituent or the like," while making full use of the advantage of "capability of achieving the 
high spatial resolution" of the conventional time gate method described above. Specifically, an object of the present 
invention is to find out a method (basic relation) for describing the behavior of light inside scattering media of different 
shapes, to find out a method for applying this basic relation to the time-resolved gate method, and thereby to provide 

45 a method and apparatus for measuring a change or an absolute value of concentration of a specific absorptive con- 
stituent inside the scattering media of various shapes by use of the foregoing relation without being affected by contours 
of scattering medium and at high spatial resolution. A further object of the present invention is to provide a measuring 
method and apparatus of absorption information inside the scattering medium, such as a tiving-body measuring method 
and apparatus, an imaging method and apparatus, a fluoroscopic method and apparatus, a tomographic image ana- 

so lyzing method and apparatus, or a mammographic method and apparatus, capable of greatly improving the measure- 
ment accuracy as to a specific absorptive constituent in a specific portion and measuring a temporal change or a spatial 
distribution thereof. 

A first measuring method of absorption information of scattering medium according to the present invention is a 
method comprising: 

55 

a first step of generating pulsed light of a predetermined wavelength; 

a second step of making said pulsed light incident in a spot shape at a light incidence position in a surface of a 
scattering medium being a measured object; 
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a third step of receiving the light having propagated inside the measured object at a photodetection position in the 
surface of said scattering medium; 

a fourth step of temporally cutting a portion out of a light signal obtained in said third step to acquire a measurement 
signal corresponding to the signal portion thus cut out; 
s a fifth step of repetitively performing said first to fourth steps at each of plural different timings and deriving a time 

integration value and a mean optical pathlength for each of plural said measurement signals obtained at said plural 
timin'^s; and 

a sixth step of calculating a change amount of absorption coefficient occurring between said plural timings, based 
on a predetermined relation among plural said time integration values, plural said mean optical pathlengths. and 
10 the change amount of absorption coefficient. 

In this case, the method may further comprise a seventh step of calculating a change amount of concentration of 
an absorptive constituent occurring between said plural timings, based on a predetermined relation among the change 
amount of absorption coefficient obtained in said sixth step, an absorption coefficient per unit concentration of the 

15 absorptive constituent, and the change amount ot concentration of the absorptive constituent. 

In said sixth step, a change amount of concentration of an absorptive constituent occurring between said plural 
timings may be calculated based on a predetermined relation among said plural time integration values, said plural 
mean optical pathlengths, an absorption coefficient per unit concentration of the absorptive constituent, and the change 
amount of concentration of the absorptive constituent. 

20 A second measuring method of absorption information of scattering medium according to the present invention is 

a method comprising: 

a first step ot generating pulsed light of a predetermined wavelength; 

a second step of making said pulsed light incident in a spot shape at a light incidence position in a surface of a 
25 scattering medium being a measured object; 

a third step of receiving the light having propagated inside the measured object at a photodetection position in the 
surface of said scattering medium; 

a fourth step of temporally cutting a portion out of a light signal obtained in said third step to acquire a measurement 
signal corresponding to the signal portion thus cut out; 
30 a fifth step of repetitively performing said first to fourth steps at each of plural different measurement positions 

while fixing a distance between said light incidence position and said photodetection position, and deriving a time 
integration value and a mean optical pathlength for each of plural said measurement signals obtained at said plural 
measurement positions; and 

a sixth step of calculating a difference of absorption coefficient occurring between said plural measurement posi- 
55 tions, based on a predetermined relation among plural said time integration values, plural said mean optical path- 

lengths, and the difference of absorption coefficient. 

In this case, the method may further comprise a seventh step of calculating a difference of concentration of an 
absorptive constituent occurring between said plural measurement positions, based on a predetermined relation among 
40 the difference of absorption coefficient obtained in said sixth step, an absorption coefficient per unit concentration of 
the absorptive constituent, and the difference of concentration of the absorptive constituent. 

In said sixth step, a difference of concentration of an absorptive constituent occurring between said plural meas- 
urement positions may be calculated based on a predetermined relation among said plural time integration values, 
said plural mean optical pathlengths, an absorption coefficient per unit concentration of the absorptive constituent, and 
45 the difference of concentration of the absorptive constituent. 

A third measuring method of absorption information of scattering medium according to the present invention is a 
method comprising: 

a first step of generating pulsed light of different wavelengths at which scattering coefficients of a measured object 
50 are equal to each other or are regarded as equal to each other; 

a second step of making said pulsed light incident in a spot shape at a light incidence position in a surface of a 
scattering medium being the measured object; 

a third step of receiving the light having propagated inside the measured object at a photodetection position in the 
surface of said scattering medium; 
55 a fourth step of temporally cutting a portion out of a light signal obtained in said third step to acquire a measurement 

signal corresponding to the signal portion thus cut out so as to keep an identical temporal relation for each of said 
pulsed light of the different wavelengths; 

a fifth step of obtaining a time integration value and a mean optical pathlength for each of plural said measurement 
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a light source for generating pulsed light of a predeternnined wavelength; 

a light incidence section for making said pulsed light incident in a spot shape at a light incidence position in a 
surface of a scattering medium being a measured object; 

a light receiving section for receiving the light having propagated inside the measured object at a photodetection 
position in the surface of said scattering medium; 

a signal detecting section for temporally cutting a portion out of a light signal obtained in said light receiving section 
to acquire a measurement signal corresponding to the signal portion thus cut out; 

a first arithmetic section for performing a measurement for acquiring said measurement signal at each of plural 
different measurement positions while fixing a distance between said light incidence position and said photode- 
tection position, and deriving a time integration value and a mean optical pathlength for each of plural said meas- 
urement signals obtained at said plural measurement positions; and 

a second arithmetic section for calculating a difference of absorption coefficient occurring between said plural 
measurement positions, based on a predetermined relation among plural said time integration values, plural said 
mean optical pathlengths, and the difference of absorption coefficient. 



In this case, said second arithmetic section may further calculate a difference of concentration of an absorptive 
constituent occurring between said plural measurement positions, based on a predetermined relation among the dif- 
ference of absorption coefficient obtained in said second arithmetic section, an absorption coefficient per unit concen- 
tration of the absorptive constituent, and the difference of concentration of the absorptive constituent. 

Also, in said second arithmetic section, a d Iff erence of concentration of an absorptive constituent occurring between 
said plural measurement positions may be calculated based on a predetermined relation among said plural time inte- 
gration values, said plural mean optical pathlengths, an absorption coefficient per unit concentration of the absorptive 
constituent, and the difference of concentration of the absorptive constituent. 

A third measuring apparatus of absorption information of scattering medium according to the present invention is 
25 an apparatus comprising: 

a light source for generating pulsed light of different wavelengths at which scattering coefficients of a measured 
object are equal to each other or are regarded as equal to each other; 

a light incidence section for making said pulsed light incident in a spot shape at a light incidence position in a 
surface of a scattering medium being the measured object; 

a light receiving section for receiving the light having propagated inside the measured object at a photodetection 
position in the surface of said scattering medium; 

a signal detecting section for temporally cutting a portion out of a light signal obtained in said light receiving section 
to acquire a measurement signal corresponding to the signal portion thus cut out so as to keep an identical temporal 
relation for each of said pulsed light of the different wavelengths; 
a first arithmetic section for obtaining a time integration value and a mean optical pathlength for each of plural said 
measurement signals; and 

a second arithmetic section for calculating a difference of absorption coefficient for said pulsed light of the different 
wavelengths, based on a predetermined relation among plural said time integration values, plural said mean optical 
"to pathlengths, and the difference of absorption coefficient. 

In this case, said second arithmetic section may further calculate a concentration of an absorptive constituent, 
based on a predetermined relation among the difference of absorption coefficient obtained in said second arithmetic 
section, absorption coefficients per unit concentration of the absorptive constituent for said pulsed light of the different 
-^5 wavelengths, and the concentration of the absorptive constituent. 

Also, in said second arithmetic section, a concentration of an absorptive constituent may be calculated based on 
a predetermined relation among plural said time integration values, plural said mean optical pathlengths, absorption 
coefficients per unit concentration of the absorptive constituent for said pulsed light of the different wavelengths, and 
the concentration of the absorptive constituent 

Further, the apparatus may be arranged preferably in such a way 

that said light receiving section receives the light at a plurality of photodetection positions; 
that said signal detecting section acquires measurement signals respectively corresponding to said pulsed light 
of the different wavelengths at each of said plurality of photodetection positions; 

that said first arithmetic section derives the time integration value and the mean optical pathlength for each of said 
plural measurement signals; and 

that said second arithmetic section calculates said concentration of the absorptive constituent, based on a prede- 
termined relation among plural said time integration values, plural said mean optical pathlengths. absorption co- 
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efficients per unit concentration of the absorptive constituent tor said pulsed light of the different wavelengths, and 
the concentration of the absorptive constituent. 

Also, the apparatus may further comprise a third arithmetic section for moving the light incidence position in said 
5 light incidence section and the photodetection position in said light receiving section along a periphery of the scattering 
medium being the measured object, calculating the concentration of the absorptive constituent obtained in said second 
arithmetic section for each of plural measurement positions achieved by said moving, and thereby calculating a tom- 
ographic image o1 concentration distribution in the scsttering medium. 

The present invention will be more fully understood from the detailed description given hereinbetow and the ac- 
70 companying drawings, which are given by way of illustration only and are not to be considered as limiting the present 
invention. 

Further scope of applicability of the present invention will become apparent from the detailed description given 
hereinafter. However, it should be understood that the detailed description and specific examples, while indicating 
preferred embodiments of the invention, are given by way of illustration only, since various changes and modifications 
IS within the spirit and scope of the invention will be apparent to those skilled in the art from this detailed description. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a schematic diagram to show the track of the photon having propagated inside the scattering medium, 
20 Fig, 2 is waveform diagrams for explaining measurement signals (cut-out signals) according to the present inven- 

tion. 

Figs. 3A and 3B are schematic diagrams, each showing a relation between the cut-out time and the measuring 
portion (a portion inside the scattering medium). 

Figs. 4A and 48 are schematic diagrams, each showing the relation between the cut-out timing and the measuring 
25 portion (a portion inside the scattering medium). 

Fig, 5 is a schematic diagram to show the configuration of the apparatus in the first embodiment according to the 
present invention. 

Fig. 6 is a graph to show the absorption spectra of hemoglobin and myoglobin. 
Fig. 7 is a schematic diagram to show another method for obtaining the measurement signal. 
30 Fig, 8 is waveform diagrams to show the relation among the output light, the gate signal, and the measurement 

signal In the case of the method shown in Fig. 7. 

Figs. 9A, 96, 9C and 9D are schematic diagrams, each showing a light incidence method to the scattering medium. 
Figs. 10A, 108 and 10C are schematic diagrams, each showing a light receiving method. 

Fig. 11 is a schematic diagram to show the configuration of the apparatus in the second embodiment according 
35 to the present invention. 

Fig. 12 is a schematic diagram to show the configuration of the apparatus in the third embodiment according to 
the present invention. 

Fig. 13 is a schematic diagram to show the configuration of the apparatus in the fourth embodiment according to 
the present invention. 

40 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 
(Principle of the Present Invention) 

^s The principle of the present invention will be described first. Herein, the problem in the photon diffusion approxi- 

mation normally used will be first clarified and then the principle of the present invention will be described. 

In general, the behavior of photon inside the scattering medium is analyzed by the photon diffusion equation based 
on the photon diffusion theory. Namely, the behavior of light inside the scattering medium is approximated by the photon 
diffusion equation, this equation Is solved to obtain a relation between an optical characteristic or an optical constant 

50 of the measured object and the output light, and absorption information or the like can be derived from this relation. 
However, the method based on the photon diffusion equation has the following problem. 

Specifically, boundary conditions must be set for solving the photon diffusion equation. However, since the bound- 
ary conditions vary greatly depending upon the shape of scattering medium, the photon diffusion equation needs to 
be solved by setting new boundary conditions for every change of the shape of scattering medium. In order to achieve 

55 accurate measurements. However, scattering media for which the boundary conditions can be set accurately to some 
extent are limited to very simple shapes such as an infinite space, a semi-infinite space, a circular cylinder having the 
Infinite length, a slab spreading infinitely and having a finite thickness, and a sphere. As a result, use of approximate 
boundary conditions is indispensable In measurements of living tissues having complicated shapes, which is a cause 
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to produce large measuring errors. This problem is also discussed, for example, in the recent literature; Albert Cerussi 
et al., "The Frequency Domain Multi-Distance Method in the Presence of Curved Boundaries," in Biomedical Optical 
Spectroscopy and Diagnostics, 1 996, Technical Digest (Optical Society of America. Washington DC, 1 996) pp. 24-26. 
Summarizing the above problem, "any measuring methods that can be systematically applied to scattering media of 
5 different shapes have not been developed yet and it is impossible for the conventional technologies to accurately 
measure the concentration of a specific absorptive constituent Inside the scattering media of different shapes system- 
atically." 

Therefore, the following discussion first discloses a relation that can be systematically applied to the scattering 
media of different shapes and then discloses a method for quantifying an absorptive constituent by applying this relation 
10 to the time-resolved gate method. In the apparatus and method of the present invention that can quantify internal 
information of the absorptive constituent or the like at improved spatial resolution, the spatial resolution can be selected 
by controlling the gate time and timing thereof. 

In non-invasive measurements of the inside of the scattering medium like the living tissues, the light is incident to 
the surface of the scattering medium, the light having propagated inside the scattering medium is received at another 
position in the surface to obtain a measurement signal, and information on the concentration or the like of the absorptive 
constituent contained inside is measured from the measurement signal. On this occasion, photons are strongly scat- 
tered by scattering constituents inside the scattering medium and optical paths thereof are bent in a zig-zag pattern. 
Fig. 1 shows an example of track of a photon where the photon incident at position P is received at position Q. A time- 
resolved waveform (impulse response) of the output light obtained after incidence of the impulse light to the scattering 
20 nnedium at a reference time t = 0 is composed of a plurality of photons having flown in various zig-zag paths and thus 
having various flight distances (optical pathlengths). However, each of the photons composing the output light at an 
arbitrary time t has a constant flight distance (optical pathlength) 1 = ct. where c is the velocity of light in the medium, 
and the Beer-Lambert law holds for each photon. Namely, the survival rate of each photon is given by exp(-c^3t) where 
M-a is the absorption coefficient. More specifically, when a light pulse consisting of many photons is incident at position 
25 p and received at position Q, many photons having passed through various optica! paths are detected at the time t, 
and the quantity of detected light being the sum thereof, i.e., the survival rate, is proportional to exp(-c|iat). In the above 
discussion, the velocity of light c is regarded as constant, because the macroscopic refractive index of living body is 
a constant value nearly equal to the refractive index of water. The above fact will be called hereinafter as microscopic 
Beer-Lambert law. 

30 From the above, the impulse response h(t) with incidence of the impulse light to the scattering medium at t = 0 

becomes a function of scattering coefficient 10.3, absorption coefficient ^ig, and time t, which is expressed as follows. 



35 



/7( 0=5(^1^. t)exp{'^^ct) (1.1) 
In h{t)=\n s{ix^, tyyx^ct (1.2) 

^Jnh(t)^-ct=-t (1.3) 

Here, &{\x^, t) is a response when the absorption coefficient is fi^ = 0 (i.e., a response in the case of only scattering 
present) and the exponent term exp(- ii^ct) is a term indicating attenuation due to the absorption coefficient p.^- The 
all functions are time-causal functions which become zero when t < 0. 

Although the above discussion concerns incidence of impulse light, it is obvious that the above relation also holds 
for incidence of any pulsed light the time width of which can be regarded as sufficiently shorter than the time response 
of medium. Accordingly, the following discussion concerns incidence of pulsed light the time width of which can be 
regarded as sufficiently short. 

It should be noted here that the above equations include no variables concerning the boundary conditions. This 
is a great difference from the solution of the photon diffusion equation including variables concerning the boundary 
conditions. In the above equations the boundary conditions are reflected through the function sd^g, t) in the impulse 
response h(t). Therefore, the above equations can be widely applied to scattering media having various non-reentrant 
surface shapes. 

At this time, the output signal Jdi^, ^g, t) is expressed as follows. 
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J{\i^, fia. t)=bs(^^. t) exp{-^^ct) (2.1) 
5 In J(|i^,^^, t)=\n b+\n s(^^, ty^i^ct (2.2) 

-^In ./hi II . n ^'Ct=-l (2.3) 

10 

Here, b is a coefficient proportional to the intensity of incident light. Time integration of the output signal J(Hs, ^Xg. 
t) indicates total integral output signal quantity l(^s, which corresponds to a measured value in CW spectroscopy. 
This integral output signal is given as follows from Eq. (2.1). 

75 

20 ^ J ( ^ = -i)c f "ts { fi^, t) exp ( -^l,cc) c/t (3.2) 

Now, let us consider that an arbitrary portion t = [t^, tg] (t^ < tg) is cut out of the output signal J{\i^, t). as shown 
in Fig. 2. The integral output signal, i.e., measurement signal ly(m, ^a) at this time, is expressed by the following 
equation. 



=bjl^^sl\i^. t) exp (-M^ct) dt (4 . 1) 



30 



-^JrfP^' exp (-n,ct) dt U .2) 

Further, the following is derived from Eqs. (4.1 ) and (4.2). 

cJ^Hsiii^, t) exp(-[i,ct) dt 

jl^s{\i^, t)exp{-\i^ct} dt 
^-cit^ 

Here, <tT> is a mean flight time of photons detected and Lt(^s. I^a) 's a mean optical pathlength thereof. This idea 
is first disclosed by the present invention. Next, integrating Eq. (5), the following is attained. 
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In Ir(v,»)^^) 

=-/o'^r(>^5'J»J^a+ln j\^s()^^,t)dtnn b {6) 



Here, the second and third terms in the right side are integration constants, which are obtained by setting = 0 
in Eq. (4.1). Equation (6) includes no variables concerning the boundary conditions, though it is a naturally expected 
result. This is greatly different from the solution of the photon diffusion equation always including the variables con- 
10 cerning the boundary conditions. In aforementioned Eq. (6), the boundary conditions are reflected in the measurement 
signal ii^) and mean optical pathlength \x^y Accordingly, Eq. (6) can be widely applied to the scattering 

media having various non-reentrant surface shapes. 

Here, the following can be understood by considering tracks of photons composing the measurement signal ob- 
tained by cutting a portion out of the output signal \x^, t). When only the initial part is cut out of the output signal 
15 J(^s, iig, t). photons of short optical pathlengths are detected as described previously, and thus, the detected photons 
are those having propagated in the narrow portion indicated by solid-line hatching (T^) in Fig. 3A and Fig. SB, which 
improves the spatial resolution. Fig. 3A shows an arrangement of transmission type measurement while Fig. 3B an 
arrangement of reflection type measurement. As the cut-out time T = tg - 1^ increases (T^-^Tg-^Tg), photons having 
propagated in successively increasing areas come to be detected as shown by dotted-oblique-line hatching (T2) and 
20 dotted-vertical-line hatching (T3) in Fig. 3A and Fig. SB, whereby the quantity of light (or quantity of signal) utilized in 
measurement increases, but the spatial resolution is lowered. 

Further, in the case of an intermediate portion being cut out of the output signal Jda^, jig, t), as the cut-out timing 
becomes later, the measured portion can be changed as shown in Fig. 4A and Fig. 4B (X^~^X2i- The measurement of 
this type can also be realized in measurements for obtaining a difference between two signals obtained by different 
time gates. Namely, it is utilized that when a later part is cut out of the output signal J(\x^, ix^, t), it includes more 
information about the peripheral portion. Therefore, in the present invention, the spatial resolution and the information- 
acquiring portion can be selected by controlling the foregoing cut-out time T = t2 - t-j and cut-out timing. 

It is noted that Fig. 3A. Fig. SB, Fig. 4A and Fig. 4B show ranges of tracks of ordinary photons forming a certain 
proportion of photons included in the detected fight (output signal). Accordingly, the detected photons also include 
30 those having propagated outside the illustrated ranges, of course. 

The following will describe methods for calculating information concerning absorption from measured values by 
use of foregoing Eq. (6). 



(Measurement of concentration change of absorptive constituent) 

Let us consider a case wherein the medium contains one type of absorptive constituent and the concentration 
thereof has changed so as to cause the absorption coefficient jig to change from \x^^ to With ordinary living bodies 
and scattering media, the scattering characteristics can be considered not to change with change in the concentration 
of absorptive constituent. This is just as if ink is put in milk. When s{\x^, X) or scattering coefficient is invariant before 
and after the change, the mean optical pathlength Lydi^, \x^) is given by a function of absorption coefficient by 
performing measurements while fixing the measuring position (the pulsed light incident position and photodetection 
position), and Eq. (6) holds for measured values before and after the change. Accordingly, the following equation is 
derived from Eq. (6), using ^ig^ and \ji^2 before and after the change. 



50 Then applying the mean value theorem, the following equation is attained from Eq. (7). 

MJlI , |I -) 

M^^s^ p^^) -(^a.-^^.r)M^.^^^.) (8) 



However, 's a suitable value satisfying < 1^ < or \l^^>[x^> ^^2. The above tells us that once the mean 
optical pathlength L-|-(Hs, p^) '® *^nown, the difference between the absorption coefficients before and after the change, 
M^a2 - 1^1. can be calculated from the values of measurement signal l-r(Us, Pa) before and after the change. 
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In this case, the mean optical pathlength Lj(\i^, p^) can be further expressed as follows using coefficient a. 

5 

However, a is a suitable value satisfying 0 < a < 1. In this case, Lj{\i^, ^) is a monotonic function and, because 
values thereof at and at normally nearly equal to each other, a = 1/2 may be assumed. The mean optical 

^i>.»iiiwt.yaiw *-*]'\^*e, ""^ '-Tv^*s' f*aZ^ ^' " * — h-^" — ^-^^ ■ ° / - a '/ 

form of measurement signal (mean time delay), as indicated in foregoing Eq. (5). 

Namely, values are calculated of the equation in the second line of aforementioned Eq. (5). Here, since c is the velocity 
of light in the medium, it can be known by a well-known method or another method. The term included in the numerator 
and the denominator, which is s(^s» t)exp{-|iaCt) = A, is the output signal Jdi^, i^a. t) in foregoing Eq. (2.1 ) multiplied by 
1/b, This coefficient b can be swept out of the integration and, because both the numerator and denominator include 
the coefficient, it is eliminated. This means that the above calculation obtains the same result by use of either A or the 
measured value Jdig, ^ig, t) of bA. Therefore, the numerator Is an integration of tJ(|Xg, iig, t) and the denominator an 
Integration o1 J(fi^, ^i^, t) over t each from t-, to tg. Calculation of these two terms and the ratio thereof can be done at 
high speed using a computer. Also, the above calculation of barycenter can be made not only by the method for directly 
calculating the numerator and the denominator in the second line of Eq. (5) as described above, but also by one of 
various methods equivalent thereto. In either case, the result of this calculation gives the mean flight time and it is 
2^ multiplied by aforementioned c to obtain the mean optical pathlength. 

For calculating concentration change AV of an absorptive constituent, the following equation derived from the Beer- 
Lambert law is used. 
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Here, e is an absorption coefficient (or extinction coefficient) per unit concentration of the absorptive constituent, 
which can be measured by a spectrophotometer. Absorption coefficients (or extinction coefficients) of many absorptive 
substances are published in various literature. Therefore, the following is obtained from Eqs. (8) to (10). 

^V=-,-^An!j^^^ (11) 

Equation (1 1 ) above clearly shows the method for measuring the concentration change AV of absorptive constituent 
in the scattering medium. 

(Application of measurement of concentration change of absorptive constituent) 



The above method can be applied to measurements of temporal changes of concentrations of plural absorptive 
constituents by using pulsed light of different wavelengths while fixing the measuring position. Further, the method can 
also be applied to measurements of temporal change in the concentration of hemoglobin in a predetermined portion, 
by repeating measurements at plural positions along a periphery of a measured object and performing a tomogram 
reconstruction arithmetic as in X-ray CT. optical CT. or the like. 

It is also possible to measure a distribution of difference from a reference value of concentration of an absorptive 
constituent inside a scattering medium, by performing the measurements while translationally moving or scanning the 
measuring position (the light incidence position and photodetection position) relative to the measured object held so 
as to keep its thickness constant, and using a measured value at an arbitrary position, as a reference value. The 
measurement of this type can be applied to photo-mammography for carrying out diagnosis of breast cancer 

In the above measurements according to the present invention, the spatial resolution is enhanced by shortening 
the gate time (cut-out time). It is also possible to control a region (or a portion) to be measured, by making the timing 
of gate variable. Specific application examples include a fluoroscope, optical CT, a clinical monitor utilized in surgery 
or cure, and so on, in addition to the photo-mammography. (Measurement of concentration of specific absorptive con- 
stituent) 

Next described is a measurement using two pulsed light beams of wavelengths and X^, i.e., the dual-wavelength 
spectrophotometry. Now, let us suppose that the absorption coefficient of the scattering medium containing one ab- 
sorptive constituent is yx^^ at wavelength and \x^2 wavelength A^. It Is also assumed that the scattering coefficients 
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l^si ^l^e medium at the wavelengths and are equal to each other or nearly equal to each other. Such 

conditions are realized readily by selecting the wavelengths used for measurement. Accordingly, in the case of the 
dual-wavelength spectroscopy at the fixed measuring position (the pulsed light incidence position and photodetection 
position), the following relations hold assuming s = Ms- 

5 

^i^sV^B)=^^i^s2^^a)^^i^s^\^a) (12-1) 
10 l-j{^sv Ma) =^T<Ms2' ^a) ^ ^T^^s' M^) (12.2) 

Then the following is attained from Eq. (6). 
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/|!;;i7«f».'l^a)4ia^lnA (13) 



Then the following equation is obtained in the same manner as in the case of Eq. (8). 

i>^ 



In 



T^f^^^a2'\^a1)^^(^^s^ Mj +ln^^ (14) 



However, is a suitable value satisfying ^Xgi < < ^^ or \i^^>\i^> When bi/b2 = 1. that is, when incident 
light intensities of the pulsed light beams of two wavelengths are equal to each other, this Eq. (14) becomes equal to 
aforementioned Eq. (7). The coefficient b/ba can be set to b.,/b2 = 1 by adjusting the intensities of incident pulsed light 
beams. Further, it is also possible to estimate the value of b^/bg by actually measuring intensities of light sources or 
pulsed light beams. Based on the above, the difference between the absorption coefficients, \x^2 ' Mai » can be calculated 
from the values of mean optical pathlength Ljdig, p^) and b^it^. and the values of measurement signals and Ij 
(7^) obtained in the measurement using the pulsed light beams of the wavelengths and X^. 

The mean optical pathlength Ly(Ps, can be expressed as follows using the same coefficient a as described 
above. 

'-t(Ms. Mx) = OcLt(Ms, Mai ) + (1 - a)MMs. Mas) (15) 

However, a is a suitable value satisfying 0 < a < 1 . Also in this case, considering the fact that Lj is the monotontc 
function, a = 1/2 may be assumed normally. Each of the mean optical pathlengths Ly(Ms, Mai) and Lt(Hs, ^^2) can be 
obtained by calculating the barycenter of the waveform of measurement signal (mean time delay) as shown in foregoing 
Eq. (5). 

The concentration V of the specific absorptive constituent is calculated from the following equation, using the 
absorption coefficients (or extinction coefficients) and per unit concentration of the specific absorptive constituent 
at the wavelengths and X2. 

V(e2-^l) = Ma2-Mai (16) 

Here, the values of G-i and G2 can be preliminarily measured by the spectrometer. Accordingly, the following equation 
is obtained. 
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Thus, the absolute concentration V of absorptive constituent can be nneasured in the same manner as the meas- 
urement of concentration change of absorptive constituent described previously. In this case, concentrations o1 (n - 1 ) 
types ot absorptive constituents can be measured using pulsed light beams of n wavelengths of two or more kinds. 

(Two-point dual- wavelength measurement) 

When the above-stated dual-wavelength spectrophotometry is carried out at two detection distances (distances 
between the light incidence and photodetection positions), the aforementioned coefficient b^/b2 can be eliminated. 
Namely, the coefficient b^/bg is eliminated by utilizing no dependence of the coefficient b^/b2 on the light incidence and 
photodetection positions. In this case, aforementioned Eq. (17) becomes as follows. 

20 Here, \t^(K^), It2(^) and Lj^ (\x^, m^), Lt2(|^s. Mx2) are measurement signals and mean optical pathlengths obtained 

at photodetection distances 1 and 2. respectively. Further, m^-, is a suitable value satisfying ^q-, < < |i.a2 iigi - 1^1 
> and 's a suitable value satisfying \x^^ ^ ^x2 ^ »^a2 ^ Mx2 ^ ^az^ each of which can be obtained in the 

same manner as in Eq. (1 5). Further, it is a matter of course that the above-stated method can be expanded to multi- 
wavelength spectrophotometry using light beams of two or more wavelengths. 
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(Measurement of spatial distribution of concentration of absorptive constituent) 



A spatial distribution of concentration ot absorptive constituent can be measured by performing the above-stated 
measurement at multiple points. In this case, the detection distance (the distance between the light incidence and 

30 photodetection positions) at each measuring point may differ, because a measured value ot concentration of absorptive 
constituent is measured independently of other values. Namely, the method of the present invention can measure 
concentrations of the absorptive constituent inside the scattering media having various non-reentrant contours at high 
spatial resolution. On this occasion, the spatial resolution is controlled by the aforementioned cut-out time T = t2 - 1-, 
and the timing thereof. Specific application examples include the photo-mammography, fluoroscope, optical CT, and 

3£ so on. These applications make use of methods of light reception at plural points, scanning of light incidence position 
and photodetection position, time-sharing measurement, and so on. The image reconstruction arithmetic as in CT is 
carried out as occasion may demand. The feature of these measurements is the capability of quantitatively measuring 
the concentration distribution of specific absorptive constituent at high spatial resolution, as described above. 

40 [ElVlBODIMENTS] 

Embodiments of the present invention will be described with reference to the accompanying drawings. It is, how- 
ever, noted that in the following description the same elements will be denoted by the same reference symbols and 
redundant description will be omitted. 



Embodiment 1 



Fig. 5 shows the first embodiment of the apparatus of the present invention for carrying out the method of the 
present invention and illustrates the configuration ot apparatus 1 for measuring a temporal change of concentration of 

so an absorptive constituent inside the scattering medium 2. This configuration of apparatus 1 is arranged to measure a 
concentration change of an absorptive constituent in scattering medium 2 containing one type of absorptive constituent 
or to measure a concentration change of an absorptive substance comprised of a plurality of constituents. In this 
configuration, pulsed light of a sufficiently narrow time width and of predetermined wavelength X is made incident at 
position P (light incidence position) in the surface of scattering medium 2 and the light having propagated inside the 

SB scattering medium 2 is received at another position Q (photodetection position) in the surface. Then this measurement 
is repeated to quantify a change of concentration of the absorptive constituent in a relatively narrow portion inside the 
scattering medium 2. In this case, the change of concentration of the absorptive constituent can be quantified by using 
a concentration of the absorptive constituent obtained in the first measurement, as a reference value. The measuring 
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apparatus 1 is integrally housed in one casing. 

Light source 10 is a laser diode or the like and generates the pulsed light of wavelength X. In this case, the wave- 
length is selected depending upon the scattering medium 2 or the absorptive constituent to be measured. In measure- 
ments of living bodies, oxygenated and deoxygenated hemoglobin and oxygenated and deoxygenated myoglobin Is 
s often measured and absorption spectra of those absorptive constituents are shown in Fig. 6. Therefore, the light of 
600 nm to 1.3 \irr\ is normally used for the measurements of living bodies. The pulse width is normally approximately 
20 ps to 200 ps. The light source may be a light-emitting diode, a pulse laser, or the like, as well as the laser diode. 

The pulsed light emitted from the light source 1 0 is made incident through light guide 1 2 to the surface of scattering 
medium 2 being a measured object. The space between the light guide 12 and the scattering medium 2 is very small 
10 in the embodiment of Fig. 5. However, in practice, this space may be expanded and may be filled with a liquid substance 
or a jelly substance (hereinafter referred to as an interface material) having the refractive index and scattering coefficient 
nearly equal to those of the scattering medium 2. Namely, no problem arises, because the light propagates in this 
interface material to enter the measured object. In the case wherein reflection on the surface of scattering medium 2 
is problematic. Influence of the surface reflection or the like can be decreased by properly selecting the interface 
material. 

The light having propagated inside the scattering medium 2 is received by light guide 13 placed at position Q 
(photod election position) the distance r apart from the aforementioned light incidence position R The interface material 
may also be used herein for the same reason as above. Photodetector 14 converts the light signal to an electric signal, 
amplifies it if necessary, and outputs it. The photodetector 14 may be selected from a phototube, a photodiode, an 

20 avalanche photodiode, a PIN photodiode, a streak camera, a streak scope, a photo-oscilloscope, and so on, as well 
as a photomultiplier tube. In selection of the photodetector 1 4, it is expected to have the spectral sensitivity character- 
istics for detecting the light of predetermined wavelengths and the necessary temporal response speed. For weak light 
signals, a high-gain photodetector is used. Further, the time correlation photon counting method for counting photons 
may be applied. The places other than the light receiving surface of the photodetector are desirably constructed in 

25 structure for absorbing or intercepting the light. 

Signal detecting section 16 cuts a predetermined time zone out of the detection signal to output a measurement 
signal. Specifically, a gate circuit well known is utilized. On this occasion, a gate signal for cutting out the predetermined 
time zone is generated in signal processing section 11. In this case, the signal processing section 11 uses a signal 
synchronized with the pulsed light generated from the light source 10 as occasion may demand. The time resolution 

30 can be controlled by making the time width (gate time) of the gate signal variable. It is also possible to select a portion 
to be measured by changing the timing of gate signal. For example, when a fore portion of the detection signal is 
extracted by a narrow gate, the light having propagated through a narrow portion 2a inside the scattering medium is 
measured, while the light having propagated through a wider portion 2b than the portion 2a is measured with increasing 
gate time. Further, it is also possible to employ the controls as shown in foregoing Figs. 3A and 3B and Figs. 4A and 4B. 

35 ,The above acquisition of measurement signal may be made by another method, lor example, by direct gating by 

the photodetector, as shown in Fig. 7 and Fig. 8. For example, in the case of use of a photomultiplier tube, the gate 
operation can be made by applying a pulsed voltage to the photocathode, to the first dynode or another dynode, or to 
the cathode. With an avalanche photodiode, it is also possible to effect the gate operation by applying the pulsed 
voltage thereto. The streak camera also has the time gate function. In these cases, one device can carry out the 

40 photoelectric conversion and gate operation. 

First arithmetic section 16 calculates the time integration value Ijand mean optical pathlength Lj (corresponding 
to the barycenter of waveform (mean delay time)) from the measurement signal. On this occasion, the mean optical 
pathlength Lj is calculated according to aforementioned Eq. (5). Then the above measurement is repetitively carried 
out at different times. In the following discussion we consider the m-th and (m + 1)-th measurements. 

-^5 Second arithmetic section 17 substitutes the two types of aforementioned time integration values \j ^ and l-|-, m+i 

obtained in the m-th and (m -i- 1)-th measurements and the mean optical pathlength Lj(\Xq, [iy) obtained from the two 
types of aforementioned mean optical pathlengths Lj^and Lj^^^ by use of Eq. (9), into aforementioned Eq. (8) to 
calculate a change amount between the absorption coefficients of the portion of scattering medium 2, ^a(m+i> " M-am 
(primary Information), and further calculates a change amount of the absorptive constituent (secondary information) 

50 by use of aforementioned Eq. (11). At this time, sufficient accuracy is assured with a = 1/2 in aforementioned Eq. (8) 
for calculation of the mean optical pathlength Lj(|ig, [ly). These arithmetic processes are carried out at high speed by 
microcomputers or the like incorporated in the first and second arithmetic sections. 

Display/recording section 18 has a function to store the concentration information of absorptive constituent ob- 
tained as described above, and displays or records these. 

55 The above description concerned use of the pulsed light of one wavelength, but in practice pulsed light of two or 

more wavelengths can also be applied. Further, it is also possible to make the pulsed light incident at one light incidence 
position and to detect the propagating light at two or more photodetection positions. These may be detected in parallel 
or in time division. 
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The means tor making the light incident to the scattering medium 2 may be selected instead of the light guide 12 
shown in Fig. 5 from a method with a condenser lens (Fig. 9A). a method using an optical fiber (Fig. 9B), a method 
utilizing a pinhole (Fig. 9C). a method for making the light incident from inside a body like a gastrocamera (Fig. 9D), 
and so on. A thick beam of light may be made incident to the scattering medium. This case may be considered as an 
s arrangement of plural spot light sources being arrayed. 

The means for detecting the light having diffuse-propagated inside the scattering medium 2 may be either one of 
a method for direct detection (Fig. 10A), a method using an optical fiber (Fig. 10B), a method using a lens (Fig. IOC), 

anO bO on, in auOlUUH lU Ult:; IIICUll^ uoil l\j iiicf iiJiJiii ywiww- witw*. a- 

10 Embodiment 2 

While the relative positional relation is kept constant (fixed) between the light incidence position P and the photo- 
detection position Q of the pulsed light in the above first embodiment, measurements can be carried out by scan along 
the scattering medium 2 and using a concentration of the absorptive constituent at an arbitrary position as a reference 

75 value, thereby measuring a spatial distribution of difference in the concentration from the reference value. In this case, 
similarly as in first embodiment, the spatial distribution of difference in the concentration of absorptive constituent from 
the reference value can be measured by use of Eqs. (8) to (11). 

Fig. 11 shows the second embodiment of the apparatus of the present invention for carrying out the method of the 
present invention and illustrates the configuration of apparatus 1 for measuring the spatial distribution of concentration 

20 of absorptive constituent inside the scattering medium 2 like the breast lightly nipped so as to have parallel and fiat 
surfaces. In Fig. 11 the same reference symbols are used for the elements having the same functions as those shown 
in Fig. 5 associated with the foregoing first embodiment. The pulsed light of predetermined wavelength X is incident to 
the surface of scattering medium 2 and the light having propagated inside the scattering medium 2 is received at a 
position in the surface on the opposite side. On this occasion, measurements are carried out while moving the light 

25 incidence position and photodetection position of the pulsed light with keeping the relative positional relation between 
them constant. Then, the spatial distribution of concentration difference of absorptive constituent can be measured, 
for example, by using a concentration of the absorptive constituent upon a measurement at the first measuring position 
(the first light incidence position and first photodetection position), as a reference value (i.e., a reference value for 
measurement of change amount). 

30 This second embodiment is provided with first mechanical section 30 for lightly nipping the scattering medium 2 

in parallel. Namely, the scattering medium 2 like the breast is measured in a slightly flattened state. This first mechanical 
section 30 is equipped with second mechanical section 31 for carrying out measurements while synchronously moving 
the light incidence position (incidence section) and photodetection position (light receiving section) of the pulsed light 
along the scattering medium with keeping the relative positional relation constant between the light incidence position 

35 and the photodetection position of the pulsed light. Then this second mechanical section 31 generates a position signal 
indicating a scanning position and this position signal is supplied to the display/recording section 18 to be used for 
calculation and display/recording of spatial distribution. Further, wavelength selector 19 is provided in the post-stage 
of the light source 1 0 for generating the pulsed light, thereby permitting selection of pulsed light of a desired wavelength 
as occasion may demand. The other portions are the same as those in the apparatus of the first embodiment. 

40 This apparatus 1 of the second embodiment is arranged to extract the initial portion of output light by the gate 

signal, and the spatial resolution can be enhanced by narrowing the gate time on that occasion. Although the above 
arrangement used the pulsed light of one wavelength, pulsed light of two or more wavelengths may be applied in 
practice. Further, it is also possible to make the light incident at one light incidence position and to detect the propagating 
light at two or more photodetection positions simultaneously or in time division. 

45 

Embodiment 3 

Fig. 12 shows the third embodiment of the apparatus of the present invention for carrying out the method of the 
present invention and illustrates the configuration of apparatus 1 for measuring a concentration distribution of absorptive 

50 constituent inside the scattering medium 2. In this case, the scattering medium 2 like the breast is nipped lightly In 
Fig. 12 the same reference symbols are used for the elements having the same functions as those shown in Fig. 5 
associated with the first embodiment and shown in Fig. 11 associated with the second embodiment. This configuration 
of apparatus 1 employs pulsed light beams of n different wavelengths to at which the scattering medium 2 shows 
scattering coefficients equal to each other or regarded as equal to each other, and can perform imaging of concentra- 

55 tions of (n - 1) types of absorptive constituents by use of aforementioned Eq. (17). This third embodiment has the 
configuration very similar to the aforementioned second embodiment, but the measuring parameter is the spatial dis- 
tribution of concentration of absorptive constituent. In addition, the measuring position (the light incidence position and 
photodetection position) is moved, but the relative positional relation between the light incidence position (light inct- 



15 

BNSDOCID: <EP 0834277A1 I > 



EP 0 834 277 A1 



dence point) and the photodetection position (ligint receiving point) of the pulsed light may change, because a concen- 
tration of the absorptive constituent at each scan position is nneasured based on foregoing Eq. (17). Specifically, the 
second ennbodinnent described above was arranged to pinch the scattering medium 2 in parallel and to scan it with 
keeping the relative positional relation constant between the light incidence position and the photodetection position, 
s whereas the third embodiment has no such limitation. The reason is as described in the foregoing section of (meas- 
urement of spatial distribution of concentration of absorptive constituent). 

Light sources 10^ to 10^ are laser diodes, which successively generate pulsed light beams of n different wave- 
lengths. On this occasion, oscillation timings of light sources 1 0-, to 1 0^ are controlled by a signal from signal processing 
section 1 1 . The pulsed light beams from the light sources 1 0^ to 1 0^ are multiplexed in multiplexer 20 and the multiplexed 

10 light is made incident trough light guide 12 to the surface of scattering medium 2 being a measured object. 

The space between the light guide 12 and the scattering medium 2 is very small in the embodiment of Fig. 12. 
However, in practice, this space may be expanded and may be filled with the interface material having the refractive 
index and scattering coefficient nearly equal to those of the scattering medium 2, as in the first embodiment. It poses 
no problem at all, because the pulsed light propagates in this interface material to enter the measured object. When 

fs reflection on the surface of scattering medium 2 is problematic, influence of the surface reflection or the like can be 
reduced by properly selecting the interface material. 

The light having propagated inside the scattering medium 2 is received by light guide 13 placed at the position on 
the opposite side to the light incidence position. The interface material may be used herein for the same reason as 
above. A light signal from the light guide 13 is subjected to wavelength selection by wavelength selector 21 and the 

20 light signal corresponding to the incident pulsed light beam of each wavelength is guided through the light guide 1 3 to 
next-stage photodetector 14. 

The scattering medium 2 is nipped lightly by the first mechanical section 30 as in the second embodiment. This 
first mechanical section 30 is provided with second mechanical section 31 for synchronously moving the light incidence 
position and photodetection position of the pulsed light. This second mechanical section 31 generates a position signal 

2S indicating a scanning position and this position signal is supplied to the display/recording section 18 to be used for 
calculation of spatial distribution and display/recording of result thereof. 

The photodetector 1 4 converts the light signal of each wavelength described above to an electric signal, amplifies 
it if necessary, and outputs a detection signal. The photodetector 14 may be selected from the phototube, photodiode, 
avalanche photodiode, PIN photodiode, streak camera, streak scope, photooscilloscope, and so on, In addition to the 

30 photomultiplier tube. In selection of the photodetector 14, it is expected to have the spectral sensitivity characteristics 
for detecting the light of predetermined wavelengths and the necessary temporal response speed. For weak light 
signals, a high-gain photodetector is used. Further, the time correlation photon counting method for counting photons 
may be applied. The other places than the light receiving surface of photodetector 14 are desirably constructed in the 
structure for absorbing or intercepting the light. 

55 Signal detecting section 15 cuts a predetermined time zone out of the aforementioned detection signal to output 

a measurement signal. A gate signal for cutting out the predetermined time zone is generated in signal processing 
section 11. The spatial resolution can be controlled by arranging the time width (gate time) of the gate signal variable 
and the spatial resolution is enhanced by narrowing the gate time. For example, when a fore portion of the detection 
signal is extracted by a narrow gate, the light having propagated in a narrow portion 2a inside the scattering medium 

40 will be measured. With increasing gate time the light having propagated in a wider portion than the portion 2a will be 
measured, which lowers the spatial resolution. The above acquisition of measurement signal may be done by another 
method, for example, by a method of direct gating by the photodetector, as shown in Fig. 7 and Fig. 8 described 
previously. 

First arithmetic section 16 calculates the time integration value l-j- and mean optical pathlength Lj (equivalent to 
45 the barycenter of waveform (mean time delay)) from the measurement signal. On this occasion, the mean optical 
pathlength Ljis calculated according to aforementioned Eq. (5). Second arithmetic section 17 calculates concentrations 
of the (n - 1) types of absorptive constituents contained in the portion 2a inside the scattering medium 2 by use of the 
n time integration values Ij^ to Ijn, n mean optical pathlengths Lj^ to Lj^, obtained in the above measurement, and n 
extinction coefficients. These arithmetic processes are carried out at high speed by microcomputers or the like incor- 
50 porated in the first and second arithmetic sections. 

Then the above sequential measurement is repeated with changing the measuring position by the second me- 
chanical section 31 , thereby calculating concentrations of the absorptive constituents at each measuring position (the 
light incidence position and photodetection position). The display/recording section 18 has a function to store concen- 
trations of the absorptive constituents obtained as described above and displays or records these. 
55 The above case was described with the method for making the pulsed light beams of the n wavelengths succes- 

sively incident, but it is also possible to employ a method for simultaneously turning the pulsed light beams of n wave- 
lengths on and making them incident coaxially. In this case there are also applicable methods including a method for 
selecting the wavelength by wavelength selector 21 provided immediately before the photodetector 14 and a method 
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for splitting the detected light into n, effecting wavelength selection thereon, and thereafter delecting them by n pho- 
todetectors in parallel. 

Fig, 12 showed the nearly parallel nipping arrangement of scattering medium 2, but the thickness thereof may 
differ in practice, as described previously. In this case, however, the second mechanical section 31 is arranged to move 
5 along the surface of scattering medium 2. Further. Fig. 12 shows the upper and lower (vertical on the plane of figure) 
arrangement of the light incidence position and photodetection position, but image reconstruction like CT may be made 
including data measured with keeping these positions in an oblique positional relation. 

Embodiment 4 

10 

Fig. 13 shows the fourth embodiment of the apparatus of the present invention for carrying out the method of the 
present invention and illustrates the configuration of apparatus for measuring a concentration distribution of absorptive 
constituent in a cross section of scattering medium 2. The light incidence position (light incidence section) and the 
photodetection position (light receiving section) are moved along a periphery of scattering medium 2. In Fig. 13 the 

IS same reference symbols are used for the elements having the same functions as those shown in Fig. 5 associated 
with the first embodiment, shown in Fig. 11 associated with the second embodiment, and shown in Fig. 12 associated 
with the third embodiment. This configuration of apparatus uses pulsed light beams of two different wavelengths 
and X2 at which the scattering coefficients are equal to each other or are regarded as equal to each other. Then a 
concentration of a specific absorptive constituent at each measuring position (the light incidence position and photo- 

20 detection position) is obtained by use of aforementioned Eq. (17) and a distribution of concentration of the specific 
absorptive constituent in the cross section of scattering medium 2 is obtained from concentration values of the specific 
absorptive constituent at a plurality of measuring positions. This apparatus of the fourth embodiment has the configu- 
ration very similar in the other portions than around the light incidence position and photodetection position to the third 
embodiment described above, but a measurement result is obtained as a tomographic image. 

25 The scattering medium 2 being the measured object is set in circular rotation mechanical section 32 surrounding 

the scattering medium 2. Interface material 33 fills the gap between the scattering medium 2 and the internal wall of 
rotation mechanical section 32 and the rotation mechanical section 32 is arranged to rotate while fixing the scattering 
medium 2, The rotation mechanical section 32 is equipped with optical fiber 12 for incidence of light and light guide 13 
for reception of light. This rotation mechanical section 32 generates a position signal indicating a scanning position 

30 and this position signal is supplied to third arithmetic section 22 to be used for reconstruction of tomographic image. 

The. light source 10 is a laser diode or the like to generate two pulsed light beams of different wavelengths and 
wavelength switch 1 9 lets the two pulsed light beams pass alternately as switching them one from another. Oscillation 
of the light source 10 and switch timing of wavelength switch 19 are controlled by a signal from signal processing 
section 1 1 . On this occasion, it is also permissible to employ pulsed light sources of two different wavelengths and turn 

35 them alternately on. It is also possible to employ a method for making the pulsed light beams of two wavelengths 
simultaneously incident. In this case, the wavelength selector 19 is omitted and a wavelength selector is provided in 
a preceding stage of photodetector 14 described below to select the wavelength used for measurement. 

The pulsed light emitted from the wavelength switch 19 is incident through the optical fiber 12 to the surface of 
interface material 33 surrounding the scattering medium 2 being the measured object. The interface material 33 is the 

40 one already described in the first embodiment. The light having propagated inside the scattering medium 2 and interface 
material 33 is received by the light guide 1 3 disposed at the position on the opposite side to the light incidence position. 
A light signal from the light guide 13 is guided to the next -stage photodetector 14. 

The photodetector 14 converts light signals of the two wavelengths to electric signals, amplifies them if necessary 
and outputs detection signals. The photodetector 14 may be selected from the phototube, photodiode, avalanche 

45 photodiode. PIN photodiode, streak camera, streak scope, photo-oscilloscope, and so on, in addition to the photom- 
ultiplier tube. For weak light signals, a high-gain photodetector is used. Further, the time correlation photon counting 
method for counting photons may be applied. The other pjaces than the light receiving surface of photodetector 14 are 
desirably constructed in the structure for absorbing or intercepting the light. 

Signal detecting section 15 cuts a predetermined time zone out of the detection signal to output a measurement 

50 signal. A gate signal for cutting out the predetermined time zone is generated in the signal processing section 11 . The 
spatial resolution can be controlled by arranging the time width (gate time) of the gate signal variable and the spatial 
resolution can be enhanced by narrowing the gate time. The above acquisition of measurement signal may be replaced 
by another method, for example, by a method for direct gating by the photodetector, as shown in Fig. 7 and Fig. 8 
described above. 

55 First arithmetic section 16 calculates the time integration value IT and mean optical pathlength Lj (equivalent to 

the barycenter of waveform (mean time delay)) from the measurement signal. On this occasion, the mean optical 
pathlength is calculated according to aforementioned Eq. (5). Second arithmetic section 17 calculates the concentration 
of the specific absorptive constituent at the aforementioned measuring position by use of the two time integration values 
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Itt and l-pg, the two mean optical pathlengths L-^ and Lx2, obtained in the above measurement, and the two extinction 
coefficients. 

These arithmetic processes are carried out at high speed by microcomputers or the like Incorporated in the first 
and second arithmetic sections. 

Next, the above sequential measurement is repeated with changing the measuring position (the light incidence 
position and photodetection position) by rotating the rotation mechanical section 32 by a suitable angle, thereby cal- 
culating concentrations of the absorptive constituent at respective measuring positions. The third arithmetic section 
22 calculates a concentration distribution (tomographic image) of the absorptive constituent in the cross section from 
the values of concentrations of the absorptive constituent at the plural measuring positions, obtained in this way. 

Display/recording section 1 8 has a function to store the concentration distribution of absorptive constituent obtained 
as described above and displays or records these as a tomographic image. Techniques or hardware well known in the 
field of CT can be utilized for the above reconstruction of tomographic image. 

The above embodiment was arranged to keep the relative positional relation constant between the light incidence 
position and the photodetection position of the pulsed light, but this relative positional relation may be arranged as 
variable. Since this case increases the sampling density for the scattering medium and also increases the number of 
measurement data, the spatial resolution of tomographic image is enhanced. Further, it is also possible to realize CT 
for three-dimensional analysis by three-dimensionally variably controlling the relative positional relation between the 
light incidence position and the photodetection position of the pulsed light. 

As described above, the absorption information measuring methods and apparatus of scattering medium according 
to the present invention can measure at high spatial resolution the concentration change or the absolute concentration 
of a specific absorptive constituent inside the scattering medium of an arbitrary shape comprised of non-reentrant 
surfaces. Accordingly, the present invention can greatly improve the measurement accuracy as to the specific absorp- 
tive constituent in a specific portion and can provide the measuring methods and apparatus of absorption information 
inside the scattering medium, such as the living-body measuring method and apparatus, imaging method and appa- 
ratus, fluoroscopic method and apparatus, tomographic image analyzing method and apparatus, or mammographic 
method and apparatus, capable of measuring the temporal change or the spatial distribution without being affected by 
the contours of scattering medium and at high spatial resolution. 

From the invention thus described, it will be obvious that the invention may be varied in many ways. Such variations 
are not to be regarded as a departure from the spirit and scope of the invention, and all such modifications as would 
be obvious to one skilled in the art are intended for inclusion within the scope of the following claims. 

The basic Japanese Application No. 263047/1996 filed on October 3, 1996 is hereby incorporated by reference. 
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Claims 

1. A measuring method of absorption information of scattering medium comprising: 



a first step of generating pulsed light of a predetermined wavelength; 

a second step of making said pulsed light incident in a spot shape at a light incidence position in a surface of 
a scattering medium being a measured object; 

a third step of receiving the light having propagated inside the measured object at a photodetection position 
in the surface of said scattering medium; 

a fourth step of temporally cutting a portion out of a light signal obtained in said third step to acquire a meas- 
urement signal corresponding to the signal portion thus cut out; 

a f'^th step of repetitively performing said first to fourth steps at each of plural different timings and deriving a 
time integration value and a mean optical pathlength for each of plural said measurement signals obtained at 
said plural timings; and 

a sixth step of calculating a change amount of absorption coefficient occurring between said plural timings, 
based on a predetermined relation among plural said time integration values, plural said mean optical path- 
lengths, and the change amount of absorption coefficient. 

2. A measuring method of absorption information of scattering medium according to Claim 1, further comprising a 
seventh step of calculating a change amount of concentration of an absorptive constituent occurring between said 
plural timings, based on a predetermined relation among the change amount of absorption coefficient obtained in 
said sixth step, an absorption coefficient per unit concentration of the absorptive constituent, and the change 
amount of concentration of the absorptive constituent. 

3. A measuring method of absorption information of scattering medium according to Claim 1 , wherein in said sixth 
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step, a change amount of concentration of an absorptive constituent occurring between said plural timings is cal- 
culated based on a predetermined relation among said plural time integration values, said plural mean optical 
pathlengths, an absorption coefficient per unit concentration of the absorptive constituent, and the change amount 
of concentration of the absorptive constituent. 

5 

4. A measuring method of absorption information of scattering medium comprising: 

a first step of generating pulsed light of a predetermined wavelength; 

a second step of making said pulsed light incident in a spot shape at a light incidence position in a surface of 
10 a scattering medium being a measured object; 

a third step of receiving the light having propagated inside the measured object at a photodetection position 
in the surface of said scattering medium; 

a fourth step of temporally cutting a portion out of a light signal obtained in said third step to acquire a meas- 
urement signal corresponding to the signal portion thus cut out; 
75 a fifth step of repetitively performing said first to fourth steps at each of plural different measurement positions 

while fixing a distance between said light incidence position and said photodetection position, and deriving a 
time integration value and a mean optical pathlength for each of plural said measurement signals obtained at 
said plural measurement positions; and 

a sixth step of calculating a difference of absorption coefficient occurring between said plural measurement 
20 positions, based on a predetermined relation among plural said time integration values, plural said mean optical 

pathlengths, and the difference of absorption coefficient. 

5. A measuring method of absorption information of scattering medium according to Claim 4, further comprising a 
seventh step of calculating a difference of concentration of an absorptive constituent occurring between said plural 

2S measurement positions, based on a predetermined relation among the difference of absorption coefficient obtained 

in said sixth step, an absorption coefficient per unit concentration of the absorptive constituent, and the difference 
of concentration of the absorptive constituent. 

6. A measuring method of absorption information of scattering medium according to Claim 4, wherein in said sixth 
30 step, a difference of concentration of an absorptive constituent occurring between said plural measurement posi- 
tions is calculated based on a predetermined relation among said plural time integration values, said plural mean 
optical pathlengths. an absorption coefficient per unit concentration of the absorptive constituent, and the difference 
of concentration of the absorptive constituent. 

35 7. A measuring method of absorption information of scattering medium comprising: 

a first step of generating pulsed light of different wavelengths at which scattering coefficients of a measured 
object are equal to each other or are regarded as equal to each other; 

a second step of making said pulsed light incident in a spot shape at a light incidence position in a surface of 
40 a scattering medium being the measured object; 

a third step of receiving the light having propagated inside the measured object at a photodetection position 
in the surface of said scattering medium; 

a fourth step of temporally cutting a portion out of a light signal obtained in said third step to acquire a meas- 
urement signal corresponding to the signal portion thus cut out so as to keep an identical temporal relation for 
45 each of said pulsed light of the different wavelengths; 

a fifth step of obtaining a time integration value and a mean optical pathlength for each of plural said meas- 
urement signals; and 

a sixth step of calculating a difference of absorption coefficient for said pulsed light of the different wavelengths, 
based on a predetermined relation among plural said time integration values, plural said mean optical path- 
50 lengths, and the difference of absorption coefficient. 

8. A measuring method of absorption information of scattering medium according to Claim 7, further comprising a 
seventh step of calculating a concentration of an absorptive constituent, based on a predetermined relation among 
the difference of absorption coefficient obtained in said sixth step, absorption coefficients per unit concentration 

55 of the absorptive constituent for said pulsed light of the different wavelengths, and the concentration of the ab- 

sorptive constituent. 

9. A measuring method of absorption information of scattering medium according to Claim 7, wherein in said sixth 
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step, a concentration of an absorptive constituent is calculated based on a predetermined relation among plural 
said time integration values, plural said mean optical pathlengths, absorption coefficients per unit concentration 
of the absorptive constituent for said pulsed light of the different wavelengths, and the concentration of the ab- 
sorptive constituent. 

5 

10. A measuring method of absorption information of scattering medium according to Claim 9, 

lA/KArAirt ir^ ^^^i^ 4l^tr-<^ i^^^n «Um li^U^ in ^^saa^.a^J n« « » 1. . li*. . _l * 1—4. :a : . ^ _ , 

««ttw.iwhi iit iiitiu otcp ulc; M^iii lo icL^civcu cai a piurciiiiy ui pi luiuutiitiutlun pu&9luons| 

wherein in said fourth step measurement signals are acquired respectively corresponding to said pulsed light 
of the different wavelengths at each of said plurality of photodetection positions; 

wherein in said fifth step the temporal integration value and the mean optical pathlength are derived for each 
of said plural measurement signals; and 

wherein in said sixth step said concentration of the absorptive constituent is calculated based on a predeter- 
mined relation among plural said time integration values, plural said mean optical pathlengths, absorption 
coefficients per unit concentration of the absorptive constituent for said pulsed light of the different wave- 
lengths, and the concentration of the absorptive constituent. 

11. A measuring method of absorption information of scattering medium according to Claim 8, further comprising an 
eighth step of moving the light incidence position in said second step and the photodetection position in said third 

?o step along a periphery of the scattering medium being the measured object, calculating the concentration of the 

absorptive constituent obtained in said sixth step for each of plural measurement positions achieved by said mov- 
ing, and thereby calculating a tomographic image of concentration distribution in the scattering medium. 



12. A measuring apparatus of absorption information of scattering medium comprising: 

25 

a light source tor generating pulsed light of a predetermined wavelength; 

a light incidence section for making said pulsed light incident in a spot shape at a light incidence position in a 
surface of a scattering medium being a measured object; 

a light receiving section for receiving the light having propagated inside the measured object at a photodetec- 
tion position in the surface of said scattering medium; 

a signal detecting section for temporally cutting a portion out of a light signal obtained in said light receiving 
section to acquire a measurement signal corresponding to the signal portion thus cut out; 
a first arithmetic section for performing a measurement for acquiring said measurement signal at each of plural 
different timings and deriving a time integration value and a mean optical pathlength for each of plural said 
^5 measurement signals obtained at said plural timings; and 

a second arithmetic section for calculating a change amount of absorption coefficient occurring between said 
plural timings, based on a predetermined relation among plural said time integration values, plural said mean 
optical pathlengths, and the change amount of absorption coefficient. 



40 1 3. A measuring apparatus of absorption information of scattering medium according to Claim 1 2, wherein said second 
arithmetic section further calculates a change amount of concentration of an absorptive constituent occurring be- 
tween said plural timings, based on a predetermined relation among the change amount of absorption coefficient 
obtained in said second arithmetic section, an absorption coefficient per unit concentration of the absorptive con- 
stituent, and the change amount of concentration of the absorptive constituent. 

45 

14. A measuring apparatus of absorption information of scattering medium according to Claim 12, wherein in said 
second arithmetic section, a change amount of concentration of an absorptive constituent occurring between said 
plural timings is calculated based on a predetermined relation among said plural time integration values, said plural 
mean optical pathlengths, an absorption coefficient per unit concentration of the absorptive constituent, and the 
50 change amount of concentration of the absorptive constituent. 



15. A measuring apparatus of absorption information of scattering medium comprising: 



a light source for generating pulsed light of a predetermined wavelength; 

a light incidence section for making said pulsed light incident in a spot shape at a light incidence position in a 
surface of a scattering medium being a measured object; 

a light receiving section for receiving the light having propagated inside the measured object at a photodetec- 
tion position in the surface of said scattering medium; 
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a signal detecting section for temporally cutting a portion out of a light signal obtained in said light receiving 
section to acquire a measurement signal corresponding to the signal portion thus cut out; 
a first arithmetic section for performing a measurement for acquiring said measurement signal at each of plural 
different measurement positions while fixing a distance between said light incidence position and said photo- 
s detection position, and deriving a time integration value and a mean optical pathlength for each of plural said 

measurement signals obtained at said plural measurement positions; and 

a second arithmetic section for calculating a difference of absorption coefficient occurring between said plural 
measurement positions, based on a predetermined relation among plural said time integration values, plural 
said mean optical pathlengths. and the difference of absorption coefficient. 

10 

16. A measuring apparatus of absorption information of scattering medium according to Claim 15, wherein said second 
arithmetic section further calculates a difference of concentration of an absorptive constituent occurring between 
said plural measurement positions, based on a predetermined relation among the difference of absorption coeffi- 
cient obtained in said second arithmetic section, an absorption coefficient per unit concentration of the absorptive 
IS constituent, and the difference of concentration of the absorptive constituent. 



17. A measuring apparatus of absorption information of scattering medium according to Claim 15, wherein in said 
second arithmetic section, a difference of concentration of an absorptive constituent occurring between said plural 
measurement positions is calculated based on a predetermined relation among said plural time integration values, 
^0 said plural mean optical pathlengths, an absorption coefficient per unit concentration of the absorptive constituent, 

and the difference of concentration of the absorptive constituent. 



18. A measuring apparatus of absorption information of scattering medium comprising: 

25 a light source for generating pulsed light of different wavelengths at which scattering coefficients of a measured 

object are equal to each other or are regarded as equal to each other; 

a light incidence section for making said pulsed light incident in a spot shape at a light incidence position in a 
surface of a scattering medium being the measured object; 

a light receiving section for receiving the light having propagated inside the measured object at a photodetec- 
30 tion position in the surface of said scattering medium; 

a signal detecting section for temporally cutting a portion out of a light signal obtained in said light receiving 
section to acquire a measurement signal corresponding to the signal portion thus cut out so as to keep an 
identical temporal relation for each of said pulsed light of the different wavelengths; 

a first arithmetic section for obtaining a time integration value and a mean optical pathlength for each of plural 
35 said measurement signals; and 

a second arithmetic section for calculating a difference of absorption coefficient for said pulsed light of the 
different wavelengths, based on a predetermined relation among plural said time integration values, plural 
said mean optical pathlengths, and the difference of absorption coefficient. 



40 19. A measuring apparatus of absorption information of scattering medium according to Claim 18, wherein said second 
arithmetic section further calculates a concentration of an absorptive constituent, based on a predetermined rela- 
tion among the difference of absorption coefficient obtained in said second arithmetic section, absorption coeffi- 
cients per unit concentration of the absorptive constituent for said pulsed light of the different wavelengths, and 
the concentration of the absorptive constituent. 

45 

20. A measuring apparatus of absorption information of scattering medium according to Claim 18, wherein in said 
second arithmetic section, a concentration of an absorptive constituent is calculated based on a predetermined 
relation among plural said time integration values, plural said mean optical pathlengths, absorption coefficients 
. per unit concentration of the absorptive constituent for said pulsed light of the different wavelengths, and the con- 
50 centration of the absorptive constituent. 



21. A measuring apparatus of absorption information of scattering medium according to Claim 20, 

wherein said light receiving section receives the light at a plurality of photodetection positions; 
55 wherein said signal detecting section acquires measurement signals respectively corresponding to said pulsed 

light of the different wavelengths at each of said plurality of photodetection positions; 

wherein said first arithmetic section derives the time integration value and the mean optical pathlength for 
each of said plural measurement signals; and 



21 

BNSDOCID: <EP 0834277A1J_> 



EP 0 834 277 A1 



wherein said second aritlimetic section calculates said concentration of the absorptive constituent, based on 
a predeternnined relation among plural said time integration values, plural said mean optical pathlengths, ab- 
sorption coefficients per unit concentration of the absorptive constituent for said pulsed light of the different 
wavelengths, and the concentration of the absorptive constituent. 

5 

22. A measuring apparatus of absorption information of scattering medium according to Claim 19, further comprising 
a third arithmetic section for moving the light incidence position in said light incidence section and the photode- 
iection posiiion in said light receiving section along a periphery of the scattering medium being the measured 
object, calculating the concentration of the absorptive constituent obtained in said second arithmetic section for 
each of plural measurement positions achieved by said moving, and thereby calculating a tomographic image of 
concentration distribution in the scattering medium. 

23. A system for determining absorption characteristics of an object, the system comprising: 

illuminating a sample with plural Incident light signals and detecting in respect of light signals that travel through 
the sample a detected light signal, removing a portion of each detected signal to obtain a plurality of meas- 
urement signals and processing said measurement signals to obtain absorption characteristics of the object. 

20 
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Fig. 6 
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Fig. 8 
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Fig. 10A 
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Fig. 11 
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Fig. 12 
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